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We have investigated the electrical resistivity of the heavy fermion superconductor PrOs4Sbi2 
in the mixed state. We found unusual double minima in the flux-flow resistivity as a function 
of magnetic field below the upper critical field for the first time, indicating double peaks in the 
pinning force density (Fp). Estimated Fp at the peak exhibits apparent dependence on applied 
field direction; composed of two-fold and four-fold symmetries mimicking the reported angular 
dependence of thermal conductivity (k). The result is discussed in correlation with the double step 
superconducting (SC) transition in the specific heat and the multiple SC-phases inferred from the 
angular dependence of k. 



PrOs4St>i2, a member of the filled skutterudite family 
with the tetrahedral (Th) point group symmetry^ is the 
first Pr-based heavy fermion (HF) superconductor with 
a SC-transition temperature Tq = 1.85 K and an upper 
critical field Hq2 — 2.45 T. 2 Highly enhanced effective 
mass has been judged by the large specific heat jump 
AC/T C = 500 mJ/mol K 2 ^ and more directly by the 
de Haas- van Alphen experiment^ 

The Sb-NQR measurement has revealed that the tem- 
perature T dependence of nuclear-spin-lattice-relaxation 
rate shows no coherence peak just below Tq but an ex- 
ponential T dependence^ which is explainable as neither 
the conventional s-wave type nor any unconventional ones 
with the line-node gap. In the zero field specific hea.l£&l 
and the dc magnetization^ double step SC-transitions 
at T\ and T2 (< 2i) have been observed, which reminds 
us the two SC-phases in UPt3i£ Measson et al. have 
carefully determined two phase boundaries T\(H) and 
T2(H) by the ac specific heat measurement Izawa et 
al. reported angular dependence of k with two differ- 
ent symmetries across a boundary in the H — T phase 
diagram. They ascribed the symmetry change to the 
distinct SC-phases with the gap function A(k) having 
two point nodes in low fields and six point nodes in high 
fields^ Huxley et al. found flux-line distortion consis- 
tent with the two-fold symmetry in k by the small an- 
gle neutron scattering experiment^ Consistently, Chia 
et al. reported Independence in the penetration depth 
indicating the Independent superfluid density below 
0.3Tcr* from which point nodes on the Fermi surface 
have been inferred. In a recent muon spin relaxation ex- 
periment, Aoki et al. found spontaneous appearance of 
an internal magnetic field below As intermetallic 

compounds, PrOs4Sbi2 is the first superconductor with 
a time-reversal symmetry (TRS) breaking, since such 
SC-state with broken TRS has been confirmed only in 
Sr2Ru04 until now^ In discussing the origin of this at- 
tractive SC-state, the existence of field induced ordered 
phase (FIOP) close to the SC-state in the H — T phase 
diagram is inferred to be important ^ZAifiiii Recently, 
the FIOP has been ascribed to an antiferro-quadrupolar 



ordering based on neutron diffraction experiments^ 

In order to have deeper insight into such an unconven- 
tional SC-state in PrOs4Sbi2, we have investigated the 
flux-flow resistivity in the mixed state that is sensitive to 
the change in vortex state and energy dissipation due to 
quasiparticles trapped within voltex coieswS^ 

Single crystals of PrOs4Sbi2 were grown by Sb-self-fiux 
method, using high-purity elements, 4N (99.99% pure)- 
Pr, 4N-Os, and 6N-Sbii^ On the samples grown in the 
same manner and from the same raw materials, de Haas- 
van Alphen oscillations had been clearly detected^ which 
indirectly proves high quality of the samples. The electri- 
cal resistivity (p) under magnetic fields was measured by 
the conventional dc four-probe method. Error in deter- 
mining angle between the current direction and the axis 
of field rotation is a few degrees. 

Figure ^ shows the magnetic field H dependence of p 
at 1.3 K for selected values of current density j. We plot 
only the field increasing curves in this paper for clarity, 
though slight hysteretic behaviors have been found. For 
j smaller than 3.8 x 10 4 A/m 2 , the absence of flux flow 
manifests in zero resistivity up to Hq%. Above j = 3.8 x 
10 4 A/m 2 , p starts to increase at a field H m (e.g., H m = 
0.58 T for j = 3.8 x 10 4 A/m 2 , note the notation in the 
inset for j = 4.7 x 10 4 A/m 2 ) where the Magnus force 
defined as F m = \j x H m \ overcomes the pinning force. 
After reaching a peak at Hi, p drops to zero at around 
Ha = 0.70 T, suggesting an enhanced pinning force at 
around Ha- 

For larger j(= 4.7 x 10 4 A/m 2 ), the field H m , = 0.29 T 
where p starts to increase, is much smaller than that for 
j = 3.8 x 10 4 A/m 2 . After reaching a peak at il2 = 
0.39 T, p drops to zero around H B = 0.48 T. On further 
increasing H, p shows a second peak at Hi = 0.69 T, 
drops again around Ha = 0.72 T, and finally reaches the 
normal state resistivity at ilc2 = 0.96 T. Double minima 
at around Ha and Hb in p(H) suggest the successive en- 
hancement of pinning force related with so-called peak 
effect in the magnetization. A minimum in p(H) has 
been observed in many superconductors including ones 
containing rare earth ions^i In those experiments, how- 
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FIG. 1: Magnetic field dependence of electrical resistivity (p) 
in PrOs4Sbi2 for selected values of current densities j(10 4 
A/m 2 ). An expanded view for j = 4.7 is shown in the inset 
where typical field names used in the text are defined. 



ever, only a single minimum has been reported just be- 
low Hq2 until now, and such double minima structure as 
in Fig. n has not been reported before except our pre- 
liminary report to the best of our knowledge^ In or- 
der to confirm the reproducibility, we have checked the 
sample dependence of p(H) as shown in Fig. |3 where 
both the resistivity and magnetic field are normalized 
at Hc2 for the comparison. The residual resistivity ratio 
RRR= P3ook/pi.9K is slightly different among the present 
samples; 27 for S#a, 24 for S#b and 34 for S#c. The 
double peak structure was always observed with repro- 
ducible features such as the positions of Ha and Hb for 
the same crystallographic geometry; j\\ [100] and H\\ [010] 
as shown in Fig. Ufa) and (b) except minor quantitative 
difference. Also for different crystallographic geometries, 
we found basically the same double peak structure except 
some anisotropy as shown in Fig. |2tc). The effect of T 
on flux- flow resistivity is shown in Fig. [21 where both Ha 
and Hb increase with decreasing T. Such an apparent 
T-dependence rules out the matching effect caused by T- 
independent alignment of pinning centers as an origin of 
this peak effect. In the present case, so-called synchro- 
nization effect may be the most probable origin of the 
peak effect in PrOs4Sbi2, 20 though there remains an im- 
portant question why the synchronization occurs at the 
two different magnetic field strengths. 

From the critical current density jc where finite volt- 
age starts to appear in the flux flow experiment, the pin- 
ning force density can be estimated as F p = \jc x H\. 
The field dependence of the estimated F p is shown in 
Fig.H F p is 2.6 x 10 4 N/m 3 at H A = 0.70 T and 5.0 x 10 4 
N/m 3 at Hb = 0.49 T which are comparable with those of 
CeRu 2 (~ 3xl0 5 N/m 3 ) and UPt 3 (~ lxlO 4 N/m 3 )ZL22. 
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FIG. 2: Double minima structure in the flux flow resistivity, 
(a) and (b) show sample dependence for the same crystallo- 
graphic configuration, and (c) for a different crystallographic 
orientation. 
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FIG. 3: Variation of the field dependence of flux-flow resis- 
tivity p in PrOs4Sbi2 at temperature T and current density 
j(10 4 A/m 2 ). 
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FIG. 4: Field dependence of the pinning force density F p esti- 
mated from the field dependence of p and j — V characteristics 
for j || [100] and ff||[010]. 



Taking into account the flux creep phenomenon, we 
have determined F p more accurately by measuring the 
current density (j) - voltage (V) characteristics under 
constant fields. Typical j — V curves for selected magnetic 
fields including 0.48 T (H B ) and 0.70 T (H A ) are plotted 
in Fig. |5J a), where the j — V curves under 0.8 T and 
0.9 T taken from Fig. [Hare included. Under 0.48 T, for 
example, the resistive voltage is zero up to jc = 1.1 x 10 s 
A/m 2 , and starts to increase with increasing j. F p esti- 
mated from this jc is also plotted as solid circles in Fig. 01 
The figure suggests the qualitative feature of F p (H) to be 
unchanged, although there exists minor contribution of 
the flux creep phenomenon to the p(H) measurements. 
To clarify the low current density behavior, pj ph C 2 i s 
plotted as a function of j in Fig. GJb) . Note that above 
Ha, a finite and weakly j-dependent voltage appears al- 
ready at small current densities depending on H, and 
starts to increase at around j\. Namely, the resistive 
voltage consists of two components between Ha and Hc2 ■ 

The peak F p estimated at Hb is found to depend on 
the field direction as shown in Fig. Interestingly, 
the angular dependence mimics the anisotropy found in 
the thermal conductivity n(<p) as a function of the an- 
gle <f> of the applied magnetic field around the direc- 
tion of thermal flow (particularly the curve for 0.78 T 
in Fig. |21 of Ref. 11). In the SC-state, k(4>) was de- 
composed into three terms; isotropic, two-fold and four- 
fold symmetric ones, k shows four-fold symmetry above 
a critical field H* , below which it changes into two- 
fold. The anisotropy was discussed in correlation with 
that in the SC-gap. The solid curve in Fig. shows 
AF p ((j)) ~ O.4cos(20) + 0.6cos(4</>) which roughly repro- 
duces the main features of AF p (</>). Hb/Hc2 ~ 0.5 for 
the experimental condition in Fig. is larger than the 
critical field H* / H C2 « 0.4 in Ref. 11. It should be 
noted that the approximate four-fold-symmetry of F p at 
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FIG. 5: (a) Typical j — V characteristics under selected mag- 
netic fields (Ha and Hb) at 1.3 K. (b) the resistivity normal- 
ized at Hc2 against current density transferred from Fig. 
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FIG. 6: Anglular dependence of the pinning force density F p 
at Hb- To exclude the effect of angular dependence of Hb, F p 
has been estimated at Hb for each field direction. The solid 
curve shows Ak ~ 0.4cos(2<^>) + O.6cos(40). Inset shows the 
T-dependence of F p at Hb- 
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FIG. 7: H — T phase diagram for the characteristic fields Ha 
and Hb estimated from the flux flow resistivity measurements. 
Hc2, H\ T, and H pe are from the dc magnetization;- and H t h 
by the thermal conductivity4- 



Hb mimics that of k(4>) in the high field phase in Ref. 11. 
The present finding of the anisotropy in F p can be an- 
other proof of the anisotropic SC-state in PrOs 4 Sbi2, 
and gives a new insight that the anisotropy must be re- 
lated with the voltex core state as was clearly identified 
in the neutron scattering experiment. 12 However, such a 
huge magnitude of oscillation in Fp (~ 30%) compared 
to those in k(4>)~ and in a recent specific heat exper- 
imenljSa (a few percent) might not be simply ascribed 
only to the anisotropic quasiparticlc excitations due to 
the anisotropic SC-gap, and theoretical studies are highly 
required to understand the mechanism leading to such a 
strong oscillation in Fp . The inset in Fig. shows the 
T-dependence of F p just at Hb- Within this temperature 
range, F p becomes an order of magnitude enhanced with 
decreasing T, indicating the smallness of the potential 
barrier associated with the pinning^ 



Finally, the fields Ha and Hb are plotted as a func- 
tion of T in the H — T phase diagram (Fig. 0) , along with 
the related characteristic fields determined by the vari- 
ous experimental methods reported previouslji&Ai; the 
transition points H' and T' observed in dM(H)/dH and 
dM(T)/dT, and the onset field H pc of the peak effect 
in the M(H)& ffth is the phase boundary inferred from 
the k measurement iii An important feature of this fig- 
ure is the presence of Ha almost on F2 (iFj; i.e., F p has 
a peak close to the lower transition temperature. This 
fact can be a key feature to clarify the origin of the dou- 
ble peak structure in the specific heat which has not yet 
been settled, including whether it is intrinsic or not. On 
the other hand, the points Hb are roughly on the lines 
connecting H po in the magnetization measurement in- 
cluding the anisotropy, and clearly above Hth(T). In 
the dc-resistivity measurement, it is impossible to deter- 
mine F p at lower temperatures, since the stronger pin- 
ning force requires the higher current density. It should 
be necessary to utilize surface impedance experiment to 
clarify the lower temperature voltex state of this exotic 
SC-material. 

In summary, we have found the double minima struc- 
ture in the flux flow resistivity. The higher field one 
appears close to the field of the lower temperature SC- 
transition in the specific heat, while the lower field one 
appears at a field above the boundary of the two SC- 
phases suggested by the K-measurements. The two-fold 
symmetry component reported in the angular depen- 
dence of k and neutron scattering experiment has been 
also found in the angular dependence of the pinning force 
density superposed upon the four-fold symmetry one. 
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